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Communication/Navigation  Outage  Forecasting  System  (C/NOFS) 


The  goal  of  the  Communication/Navigation  Outage  Forecasting  System  (C/NOFS)  is  to 
specify  and  forecast  the  equatorial  ionosphere,  including  spread  F  and  ionospheric  scintillation. 
Scintillation  i  s  a  r  andom  fluctuation  i  n  s  ignal  pha  se  a  nd  a  mplitude  t  hat  de  velops  w  hen  r  adio 
waves  pass  through  ionospheric  electron  density  irregularities.  Irregularities  that  cause  the  most 
intense  s  cintillation  o  ccur  n  aturally,  a  fter  s  unset,  w  ithin  a  ±20°  band  of  1  atitude  a  round  t  he 
magnetic  equator.  Although  the  statistical  occurrence  of  scintillation  is  relatively  well  known 
with  respect  to  location,  time  of  day  and  solar  cycle,  the  night-to-night  variability  is  much  more 
difficult  to  forecast. 

In  order  to  develop  the  capability  to  forecast  the  onset  of  ionospheric  irregularities  that 
trigger  s  cintillation,  th  e  C  /NOFS  satellite  w  as  s  pecffically  d  esigned  with  a  s  ix  in  strument 
payload  for  making  in  situ  and  remote  sensing  measurements  of  the  ionosphere.  T  he  C/NOFS 
satellite  w  as  s  uccessfully  launched  in  A  pril  2  008  in  to  a  lo  w  e  arth  o  rbit  a  t  1  3°  inclination. 
C/NOFS  is  the  first  satellite  dedicated  to  monitoring  and  forecasting  ionospheric  conditions.  The 
instruments  me  asure  electric  fields,  p  lasma  ch  aracteristics,  n  eutral  w  inds,  an  d  t  he  s  trength  o  f 
scintillation-producing  i  rregularities.  V  arious  m  odels  and  analysis  t  ools  use  t  he  s  pacecraft 
observations  a  nd  combine  t  hem  w  ith  ground-based  an  d  o  ther  s  atellite  d  ata  t  o  ch  aracterize 
ionospheric  irregularities  that  are  likely  to  produce  scintillation. 

The  six  instruments  on  the  C/NOFS  satellite  are  as  follows: 

CERTO  -  Coherent  Electromagnetic  Radio  Tomography 

•  Built  by  Naval  Research  Laboratory  (NRL) 

•  Consists  of  a  beacon  and  antenna  that  transmit  at  three  frequencies  to  ground  receivers 

•  The  signals  measured  by  the  ground  receivers  are  analyzed  to  determine  the  amount  of 
scintillation  along  the  path  between  the  C/NOFS  satellite  and  the  ground  receiver 

•  If  t  he  s  ignals  a  re  not  di  storted,  t  hen  t  here  i  s  no  s  cintillation.  T  he  m  ore  di  storted  t  he 
signals  are,  the  more  intense  the  scintillation  is. 

•  The  d  ata  m  easured  b  y  C  ERTO  ar  e  u  sed  t  o  validate  an  d  i  mprove  p  hysics  m  odels 
developed  by  AFRL/RV  to  forecast  scintillation 

CORISS  -  C/NOFS  Occultation  Receiver  for  Ionospheric  Sensing  and  Specification 

•  Built  by  Aerospace  Corporation 

•  Consists  of  a  s  pecially-designed  Global  Positioning  System  (GPS)  receiver  and  antenna 
for  remote  sensing  of  the  ionosphere 

•  It  measures  signals  transmitted  by  the  global  constellation  of  GPS  satellites 

•  The  s  ignals  ar  e  an  alyzed  tod  etermine  t  he  n  umber  o  f  ch  arged  p  articles  al  ong  t  he  p  ath 
between  the  GPS  satellite  and  the  C/NOFS  satellite 

•  The  number  of  charged  particles  in  disturbed  regions  of  the  ionosphere  vary  depending 
on  the  severity  of  the  scintillations. 


•  The  data  measured  by  the  CORISS  instrument  are  used  as  inputs  to  validate  and  improve 
physics  models  developed  by  AFRL/RV  to  forecast  scintillation 

IVM  -  Ion  Velocity  Meter 

•  Built  by  University  of  Texas  at  Dallas  (UTD) 

•  Consists  of  two  sensors: 

o  one  sensor  has  a  square  opening  and  measures  the  speed  and  direction  of  charged 
particles  th  at  are  mo  ving  p  erpendicular  to  th  e  C/NOFS  s  atellite’s  d  irection  o  f 
motion 

o  the  other  sensor  ha  s  a  round  ope  ning  and  m  easures  t  he  s  peed  a  nd  di  rection  of 
charged  particles  that  are  moving  along  the  satellite’s  orbit  track.  This  sensor  also 
measures  the  number  of  charged  particles  in  the  region  of  the  ionosphere  that  the 
C/NOFS  satellite  is  flying  through 

•  The  speed  and  direction  of  charged  particles  in  disturbed  regions  of  the  ionosphere  are 
important  indicators  of  whether  or  not  scintillation  will  occur 

•  The  da  ta  m  easured  b  y  t  he  IVM  i  nstrument  a  re  used  a  s  i  nputs  to  va  lidate  a  nd  i  mprove 
physics  models  developed  by  AFRL/RV  to  forecast  scintillation. 

NWM  -  Neutral  Wind  Meter 

•  Built  by  University  of  Texas  at  Dallas  (UTD) 

•  Consists  of  t  wo  s  ensors.  O  ne  s  ensor  m  easures  t  he  s  peed  a  nd  di  rection  of  t  he  g  as 
(uncharged  particles)  flowing  in  the  ionosphere  in  the  direction  of  the  satellite’s  motion 
and  the  other  sensor  measures  the  gas  flow  perpendicular  to  the  satellite’s  track. 

•  The  m  otion  of  t  he  g  as  affects  t  he  di  stribution  of  c  harged  pa  rticles  i  n  t  he  i  onosphere. 
When  the  gas  flow  produces  an  irregular  distribution  of  charged  particles  then  there  is  a 
higher  probability  that  scintillation  will  occur. 

•  The  data  measured  by  the  NWM  instrument  are  used  as  inputs  to  validate  and  improve 
physics  models  developed  by  AFRL/RV  to  forecast  scintillation. 

PLP  -  Planar  Langmuir  Probe 

•  Built  by  Air  Force  Research  Laboratory,  Space  Vehicles  Directorate  (AFRL/RV) 

•  It  measures  the  number  of  charged  particles  along  the  path  of  the  satellite 

•  As  th  e  s  atellite  is  mo  ving,  th  e  P  LP  in  strument  s  weeps  up  charged  p  articles  in  its  p  ath 
much  1  ike  a  but  terfly  net  s  tuck  out  o  f  t  he  w  indow  of  a  m  oving  car  w  ould  s  weep  up 
mosquitoes. 

•  When  the  charged  particles  hit  the  metal  plate  at  the  front  of  the  instrument  they  create  an 
electrical  current  signal.  The  signal  is  measured  and  analyzed  to  detennine  the  number  of 
charged  particles. 

•  When  the  ionosphere  is  undisturbed,  the  PLP  instrument  measures  a  smoothly  varying 
signal  along  the  orbital  track.  However,  when  the  ionosphere  is  disturbed,  the  signals 
measured  by  PLP  are  rapidly  fluctuating. 

•  The  PLP  signal  fluctuations  are  an  indicator  of  scintillation.  T  he  more  rapid  the  signal 
fluctuations,  the  more  intense  the  scintillation  is. 


•  The  PLP  data  are  used  in  models  to  specify  exactly  where  scintillations  are  at  the  time  of 
the  measurement.  This  is  referred  to  as  nowcasting.  The  data  are  also  used  to  validate 
and  improve  physics  models  developed  by  AFRL/RV  to  forecast  scintillation. 

VEFI  -  Vector  Electric  Field  Instrument 

•  Built  b  y  National  Aeronautics  a  nd  S  pace  Administration/Goddard  S  pace  F  ight  C  enter 
(NASA/GSFC) 

•  Consists  of  six  booms,  each  10  m  long,  at  the  end  of  which  are  sensors. 

•  The  six  sensors  (one  per  boom)  measure  the  force  that  exists  in  a  region  between  particles 
that  ha  ve  oppos  ite  c  harges.  T  his  qua  ntity  is  c  ailed  t  he  electric  f  ield  and  t  he  V  EFI 
instrument  measures  both  the  strength  and  direction  of  this  field. 

•  The  electric  field  in  the  ionosphere  is  important  because  it  s  ets  in  motion  the  charged 
particles.  T  hese  m  otions  c  an  c  ause  t  he  i  onosphere  t  o  be  come  di  sturbed  a  nd  uns  table, 
thereby  producing  scintillations. 

•  Changes  in  the  electric  field  over  time  are  an  indicator  of  future  scintillation. 

•  The  data  measured  by  the  VEFI  instrument  are  used  as  inputs  to  validate  and  improve 
physics  models  developed  by  AFRL/RV  to  forecast  scintillation. 


The  data  collected  by  the  C/NOFS  satellite  are  unprecedented  in  ionospheric  research.  The 
scientific  pa  pers  i  ncluded  i  n  t  his  r  eport  represent  t  he  1  arge  b  readth  an  d  d  epth  o  f  t  he  s  atellite 
measurements  and  their  analysis.  Among  the  new  findings  are  broad  plasma  density  decreases 
of  up  to  an  order  of  magnitude  for  a  few  months  around  the  June  solstice  that  are  believed  to  be 
due  toe  ooling  o  f  t  he  i  onosphere  a  nd  t  hermosphere  dur  ing  t  hese  t  imes.  A  ve  ry  s  urprising 
discovery  made  from  the  C/NOFS  measurements  is  deep  plasma  depletions  at  sunrise  that  appear 
to  be  associated  with  upward  plasma  drifts.  The  AFRL  physics-based  model  PBMOD  was  used 
extensively  tom  ake  comparisons  w  ith  s  atellite  da  ta.  T  he  m  odel  w  as  f  ound  t  o  s  uccessfully 
reproduce  de  nsity  de  pletions  obs  erved  ate  arly  morning  1  ocal  t  imes  dur  ing  f  our  c  onsecutive 
orbits  of  the  satellite.  In  addition,  the  model  was  found  to  reproduce  the  observed  latitudinal  and 
longitudinal  s  cale  s  izes  o  f  th  e  p  lasma  d  epletions.  In  a  nother  uni  que  finding,  t  he  C  /NOFS 
satellite  c  aptured  da  ta  d  uring  a  hi  gh  s  peed  s  tream  ( HSS)  e  vent  i  n  w  hich  t  he  i  onosphere  w  ent 
from  qui  escent  t  o  di  sturbed.  D  uring  t  his  e  vent,  pos  t-midnight  ir  regularities  in  eluded  strong 
equatorial  plasma  bubbles. 
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[i]  During  June  2008  broad  plasma  density  decreases 
(BPDs)  were  detected  repeatedly  by  the  Planar  Langmuir 
Probe  (PLP)  on  board  the  Communication/Navigation 
Outage  Forecasting  System  (C/NOFS)  satellite.  These 
density  minima,  not  to  be  confused  with  Equatorial 
Plasma  Bubbles  (EPBs),  occurred  within  15°  of  the 
equator,  consisted  of  reductions  in  plasma  density  up  to 
an  order  of  magnitude  and  extended  across  several  degrees 
in  azimuth  along  the  orbit.  Analysis  revealed  that  the  BPDs 
occurred  nearly  daily  from  May  through  July  2008  on 
C/NOFS,  and  that  the  widest  BPDs  were  observed  in  the 
vicinity  of  the  South  Atlantic  Anomaly  (SAA).  Similar 
BPDs  simultaneous  with  the  C/NOFS  measurements  were 
observed  by  instruments  on  the  CHAllenging  Minisatellite 
Payload  (CHAMP)  and  Defense  Meteorological  Satellite 
Program  (DMSP)  satellites.  An  examination  of  plasma 
densities  observed  by  the  DMSP  satellites  over  several  years 
revealed  that  these  phenomena  were  a  frequent  occurrence 
during  (1)  the  period  around  June  solstices;  during  (2)  solar 
minimum  years;  (3)  in  the  vicinity  of  the  SAA.  Neutral 
densities  were  examined  during  periods  when  BPDs  were 
detected,  and  at  times  there  are  simultaneous  neutral 
depletions.  One  possible  explanation  is  a  decrease  in 
temperature  of  both  ions  and  neutrals  in  the  equatorial 
region  at  these  times,  consistent  with  downwelling  in  the 
ionosphere  and  thermosphere.  Measurements  of  plasma 
temperatures  on  DMSP  support  this  hypothesis. 
Citation:  Huang,  C.  Y.,  F.  A.  Marcos,  P.  A.  Roddy,  M.  R. 
Hairston,  W.  R.  Coley,  C.  Roth,  S.  Bruinsma,  and  D.  E.  Hunton 
(2009),  Broad  plasma  decreases  in  the  equatorial  ionosphere, 
Geophys.  Res.  Lett.,  36,  L00C04,  doi:10.1029/2009GL039423. 

1.  Introduction 

[2]  Between  May  and  August  2008,  the  Planar  Langmuir 
Probe  (PLP)  on  the  C/NOFS  satellite  regularly  measured 
reduced  plasma  density  on  the  nightside  in  the  equatorial 
region  lasting  several  minutes  and  occurring  over  a  broad 
range  of  longitudes.  Density  within  the  depleted  region  was 
reduced  by  up  to  an  order  of  magnitude  below  the  density 
outside  the  Broad  Plasma  Decrease  (BPD).  The  June  2008 
solstice  occurred  during  an  extremely  quiet  solar  minimum 
interval.  The  minimum  Disturbance  Storm  Time  (Dst)  index 
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value  during  the  month  was  —40  nT  and  this  occurred  on 
15  June.  There  are  several  periods  when  the  Auroral 
Electrojet  (AE)  index  reaches  1000  nT,  but  these  do  not 
correlate  with  the  nearly  daily  occurrences  of  BPDs. 

[3]  Densities  on  the  DMSP  satellites  were  examined  for 
the  same  interval,  and  BPDs  were  noted  at  or  near  equato¬ 
rial  latitudes.  The  deepest  BPDs  were  observed  close  to  the 
SAA  which  encompasses  the  entire  azimuthal  area  from  the 
west  coast  of  South  America  to  the  center  of  southern  Africa 
and  from  the  geographic  equator  to  50°  S.  Density  deple¬ 
tions  were  observed  on  DMSP  up  to  40%  below  the  ambient 
density  outside  the  depleted  area. 

[4]  We  have  studied  these  depletions  using  a  large  array  of 
satellite-based  observations  of  the  ionosphere-thermosphere 
(IT)  system  in  order  to  detennine  their  climatology.  This 
report  briefly  summarizes  the  results  of  our  study  of  IT 
observations  and  our  interpretation  of  these  observations. 

2.  Instrumentation 

[5]  The  PLP  on  C/NOFS  measures  plasma  densities, 
electron  temperatures  and  density  fluctuations  at  rates 
ranging  from  32  to  1024  samples  per  second.  In  this  paper 
we  focus  on  electron  densities  on  the  nightside.  During  the 
period  in  this  study,  17-19  June  2008,  the  satellite  altitude 
when  depletions  were  observed  varied  from  400  to  600  km, 
with  the  bulk  of  measurements  made  between  400  and 
500  km. 

[6]  Density  measurements  on  DMSP  were  made  by  the 
Special  Sensors-Ions,  Electrons,  and  Scintillation  (SSIES) 
suite.  The  measurements  used  in  this  study  come  from 
SSIES-3  flown  on  DMSP  FI 6.  We  use  output  from  the 
Retarding  Potential  Analyzer  (RPA)  which  gives  plasma 
densities  and  temperatures  every  1  s,  over  the  range  102- 
106  cm-3  with  an  accuracy  of  10%. 

[7]  Plasma  densities  on  the  CHAMP  satellite  were 
obtained  from  the  PLP  [Cooke  et  al.,  2003;  Roth,  2004], 
which  monitors  the  spacecraft  potential,  ion  density,  and 
electron  temperature.  A  one-second  sweep  is  performed 
every  15  seconds.  Neutral  densities  were  measured  by 
the  STAR  accelerometer  [Bruinsma  et  al.,  2004]  every 
10  seconds. 

[s]  Neutral  densities  on  the  Gravity  Recovery  and  Climate 
Experiment  (GRACE)  satellite  were  measured  using  super- 
STAR  accelerometers,  similar  to  the  STAR  accelerometer  on 
CHAMP,  [Cheng  et  al,  2008]  with  a  precision  an  order  of 
magnitude  greater  than  the  CHAMP  instrument.  In  this  study, 
5 -second  averaged  data  were  used. 

3.  Observations 

[9]  An  example  of  the  observations  which  triggered  this 
investigation  is  shown  in  Figure  1.  In  Figure  la,  shown  is 
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Figure  1.  C/NOFS  and  DMSP  observations  of  a  Broad  Plasma  Depletion  (BPD)  on  19  June  2008.  (a,  top)  Plasma  density 
(black  line)  in  1 -minute  resolution,  as  well  as  1 -second  resolution  (green  line).  C/NOFS  altitude  is  indicated  by  the  blue 
line,  with  the  heavy  dotted  section  showing  when  the  satellite  is  in  darkness,  (middle)  Spectrum  of  the  high-time  resolution 
density  fluctuations,  (bottom)  The  magnetic  equator  is  shown  by  a  green  trace,  and  satellite  orbit  is  indicated  by  the  red 
trace,  with  the  heavy  dotted  section  indicating  when  the  spacecraft  is  in  darkness.  The  BPD  (purple  arrow)  starts  at  2253 
UT,  with  recovery  after  2310  UT.  (b)  In  the  DMSP  FI 6  plot  the  ion  density  is  plotted  as  a  function  of  magnetic  latitude 
during  an  equatorial  crossing  at  2236  UT.  The  BPD  (purple  arrow)  occurs  at  the  equator  as  indicated. 


the  PLP  plot  from  C/NOFS  on  19  June  2008.  In  Figure  la 
(top),  plasma  density  is  indicated  by  the  black  trace  which 
shows  1 -minute  averaged  data.  1 -second  data  are  illustrated 
by  the  green  trace.  Also  shown  is  the  satellite  altitude  in 
blue  with  the  heavy  dotted  line  indicating  when  the  space¬ 


craft  is  in  darkness.  In  Figure  la  (middle),  the  color 
spectrogram  shows  the  frequency  of  the  density  fluctua¬ 
tions.  Figure  la  (bottom)  shows  the  orbit  of  the  satellite  in 
red,  the  heavy  dotted  segment  showing  when  C/NOFS  is  in 
darkness.  The  magnetic  equator  is  plotted  in  green. 
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Figure  2.  (a)  Two  consecutive  nightside  orbits  of  GRACE  on  9  June  2005.  (top)  The  first  crossing  occurred  at  2200  UT, 
20.3  MLT,  329°E.  long,  (geog.),  altitude  of  468  km;  (bottom)  the  second  at  2333  UT,  20.2  MLT,  305°E.  long,  (geog.), 
470  km  altitude,  (b)  DMSP  F16  equatorial  crossings  on  the  same  day.  (top)  The  crossing  at  2114  UT,  20.5  MLT,  343°E. 
long,  (geog.);  (bottom)  the  plasma  density  at  2253  UT,  20.6  MLT,  320°E.  long.  (geog.).  Reduced  neutral  and  plasma 
densities  in  all  four  plots  can  be  clearly  seen,  (c)  CHAMP  ion  and  neutral  densities  at  2224  UT,  17  June  2008.  Satellite 
location  is  1.23  MLT,  329  km  altitude,  40°E  long,  (geog.)  (top)  The  plasma  density  is  shown,  as  a  function  of  latitude, 
(bottom)  CHAMP  Neutral  density. 


[10]  At  approximately  2253  UT,  at  an  altitude  of  about 
475  km,  and  at  40°  E  geog.  long.,  the  plasma  density 
decreases  sharply  as  illustrated  by  the  black  trace.  Densities 
gradually  recover,  and  are  at  their  pre-depleted  level  by 
2310  UT.  The  magnetic  latitude  of  the  plasma  decrease 
varies  from  —8°  to  —21°.  In  IB,  below,  is  shown  the  DMSP 
F16  plot  of  plasma  densities  for  the  equatorial  crossing  at 
2236  UT.  The  orbit  corresponds  to  a  single  crossing  of  the 
eveningside  equatorial  ionosphere  at  an  altitude  of  840  km, 
and  magnetic  local  time  (MLT)  of  20.2.  Note  that  while  the 
C/NOFS  and  DMSP  measurements  are  nearly  simultaneous, 
there  is  considerable  spatial  separation  between  the  two 
satellites.  The  magnetic  equator  is  placed  at  the  center  of  the 
plot,  and  latitudes  from  —50°  S  to  50°  N  are  shown  in 
intervals  of  10°.  Indicated  by  the  arrow,  there  is  a  large 
reduction  in  plasma  density  in  the  equator  at  2236  UT,  at 
318°  E  geographic  longitude  (geog.  long.),  corresponding  to 
the  SAA. 

[n]  Accurate  statistics  for  C/NOFS  cannot  be  obtained 
due  to  the  orbital  variations  in  altitude,  latitude  and  longi¬ 
tude  during  the  May-July  2008  interval.  However  in  this 
period,  barring  data  gaps,  on  all  but  8  days  BPDs  were 
observed  by  the  PLP.  On  most  days,  BPDs  were  seen  on 
multiple  orbits  if  the  satellite  was  in  the  vicinity  of  the  SAA 
on  the  nightside. 

[12]  As  can  be  seen  in  Figure  1,  the  extent  of  the  depleted 
region  is  large.  In  the  C/NOFS  plot,  the  area  exceeds  60°  in 
longitude  and  15°  in  latitude.  The  DMSP  depletion  is  more 
than  15°  in  latitude.  For  this  reason,  we  are  careful  to 
describe  these  as  broad  plasma  density  decreases  (BPDs), 


not  to  be  confused  with  EPBs  which  are  orders  of  magni¬ 
tude  smaller  in  spatial  extent.  However  there  is  a  relation 
between  BPDs  and  EPBs.  BPDs  generally  originate  over  the 
SAA,  and  are  relatively  shallow  and  unstructured.  Over 
subsequent  orbits  of  C/NOFS,  they  migrate  eastward,  and 
the  depletion  becomes  narrower  and  deeper.  When  the  BPD 
reaches  the  dawn  meridian,  they  are  considerably  narrower 
with  very  deep  structures  within  the  BPD.  For  more  detail, 
we  refer  the  reader  to  de  la  Beaujardiere  et  al.  [2009], 

[13]  Having  confirmed  the  presence  of  BPDs  on  DMSP 
in  June  2008,  we  examined  its  extensive  database  for  other 
examples.  We  found  that  BPDs  occur  regularly  on  the 
nightside  in  (1)  June  solstices  which  occur  (2)  during  solar 
minima,  (3)  predominantly  in  the  SAA.  During  the  period 
from  1  May-31  July  2008,  BPDs  were  observed  on  DMSP 
F16  on  87  of  the  92  days.  In  addition,  during  intervals 
centered  on  the  June  solstices  of  1995,  1996,  2005,  2006 
and  2007  we  found  BPDs  on  DMSP  repeatedly. 

[14]  Neutral  densities  were  also  examined  for  possible 
connection  with  the  BPDs.  CHAMP  and  GRACE  acceler¬ 
ometer  data  were  studied  for  the  periods  when  large  BPDs 
were  recorded.  Examples  of  neutral  density  decreases  ob¬ 
served  on  GRACE  are  shown  in  Figure  2a.  On  9  June  2005, 
two  successive  nightside  orbits  are  shown.  In  Figure  2a 
(top),  the  equatorial  crossing  occurred  at  22  UT,  at  a 
location  of  20.3  MLT,  altitude  of  468  km,  329°  E  geog. 
long.  In  Figure  2a  (bottom),  the  equatorial  crossing  occurred 
at  2333  UT,  at  20.2  MLT  and  470  km  altitude,  305°E  geog. 
long.  In  both  orbits,  reductions  in  neutral  density  up  to  45% 
can  be  seen  around  the  equator.  The  other  12  equatorial 
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Figure  3.  (top)  Plasma  density  and  (bottom)  temperature 
from  DMSP  F16  on  9  June  2005,  with  (middle)  GRACE 
neutral  density.  Plasma  and  neutral  densities  are  both 
reduced  in  the  equatorial  region.  Lower  temperatures  at  the 
equator  clearly  indicate  cooling  at  this  time. 

crossings  on  this  day  which  do  not  occur  close  to  the  SAA 
do  not  display  this  phenomenon. 

[is]  In  Figure  2b  are  shown  two  equatorial  crossings 
made  by  DMSP  F16  on  the  same  day.  Figure  2b  (top)  shows 
measurements  made  at  21 14  UT,  at  a  343°E  geog.  long.,  and 
20.5  MLT.  Figure  2b  (bottom)  shows  the  subsequent  cross¬ 
ing  at  2253  UT,  320°E  geog.  long.,  and  20.6  MLT.  These 
crossings  are  nearly  coincident  both  in  time  and  space  with 
the  GRACE  equatorial  orbits  on  the  left  side  of  Figure  2. 
The  density  reductions  inside  the  BPDs  on  DMSP  are 
approximately  40%  of  ambient  values. 

[ 1 6]  On  17  June  2008,  on  three  consecutive  equatorial 
crossings  between  2048  and  2400  UT,  the  PLP  on  CHAMP 
detected  deep  BPDs.  The  CHAMP  orbit  crossed  C/NOFS  at 
2355  UT  which  observed  highly  structured  BPDs  at  this 
time.  During  the  2224  UT  crossing,  second  of  the  three, 
CHAMP  was  located  at  329  km  altitude,  1.23  MLT,  40°E 
geog.  long.  In  Figure  2c  are  shown  the  ion  and  neutral 
densities  during  this  crossing.  The  PLP  densities  show  a 
BPD  which  is  over  30°  across  in  latitude,  with  a  plasma 
density  that  is  more  than  an  order  of  magnitude  below 
ambient  values.  The  neutral  density  decrease  is  about  10% 
below  the  densities  outside  the  neutral  depletion.  It  is  not 
clear  whether  the  difference  between  this  example  and  that 
shown  in  Figure  2a  is  due  to  the  significantly  lower  neutral 


densities  during  2008,  or  if  this  represents  a  difference  in  the 
phenomenology. 


4.  Discussion 

[17]  The  appearance  of  BPDs  can  be  explained  in  several 
ways.  For  a  number  of  reasons  we  rule  out  the  possibility  of 
EPBs  as  the  cause.  The  climatology  of  EPBs  is  in  conflict 
with  the  climatology  of  the  BPDs  [ Gentile  et  al.,  2006]. 
During  solar  minimum,  there  are  few  EPBs  which  reach 
DMSP  altitude,  and  during  June  EPB  occurrence  in  the 
SAA  region  is  low  [ Huang  et  al.,  2002],  During  times  of 
maximum  occurrence  of  EPBs  during  the  December  sol¬ 
stices  and  equinoxes,  BPDs  have  not  been  detected.  As 
already  mentioned,  the  size  of  BPDs  exceeds  normal  EPBs 
by  a  large  factor. 

[is]  A  second  possible  cause  of  BPDs  and  neutral  deple¬ 
tions  is  a  change  in  chemistry  in  the  ionosphere  and 
thermosphere.  This  occurs  during  magnetic  storms,  when 
large  changes  in  IT  densities  occur  [Crowley  et  al.,  2006]. 
Other  systemic  variations  in  the  thermosphere  occur  as  a 
result  of  global  warming  [Akmaev  et  al.,  2006],  However 
neither  of  these  mechanisms  can  account  for  the  specific 
climatology  associated  with  the  BPDs,  nor  the  simultaneous 
decreases  in  neutral  and  ion  densities  in  the  equatorial 
region. 

[19]  Coincident  reductions  in  neutral  and  plasma  densi¬ 
ties  narrow  the  range  of  possible  causes.  Simultaneous  ion 
and  neutral  density  reductions  have  been  observed  in  EPBs 
[Bencze  et  al.,  2000]  but  as  we  have  pointed  out,  these 
BPDs  cannot  be  bubbles.  The  minima  in  ion  and  neutral 
densities  can  be  described  by  a  change  in  the  topside  scale 
height,  defined  as  H  =  kT/mg  where  k  is  Boltzmann’s 
constant,  T  is  the  temperature  of  the  species,  m  the  mass 
of  the  species  and  g  the  gravitational  constant  [Rishbeth  and 
Garriott,  1969]. 

[20]  A  change  in  plasma  and  neutral  density  can  be  due  to 
a  change  in  T,  the  plasma  or  neutral  temperature,  or  m,  the 
mass  of  the  species,  or  both.  An  examination  of  the  fraction 
of  light  ions  (H+)  during  the  DMSP  observations  shows  that 
during  the  June  2005  events,  this  fraction  is  less  than  50% 
in  the  minimum  of  the  BPD,  and  the  opposite  is  true  during 
June  2008,  when  the  solar  cycle  is  in  a  deep  minimum.  This 
rules  out  a  persistent  change  in  ion  composition. 

[21]  We  suggest  that  the  most  likely  explanation  is  a 
cooling  or  downwelling  of  the  equatorial  region  during 
these  periods.  This  would  account  for  simultaneous 
decreases  in  both  neutrals  and  plasma  density.  To  verify  this 
hypothesis,  we  examined  ion  temperatures  (T,)  on  DMSP 
during  the  periods  when  BPDs  were  seen.  In  Figure  3 
(bottom),  we  show  T;  obtained  from  DMSP  F16  during 
the  2253  UT  equatorial  crossing  on  9  June  2005,  together 
with  ion  density  (Figure  3,  top)  measured  simultaneously 
and  the  GRACE  neutral  density  (Figure  3,  middle)  from  the 
2333  UT  crossing.  There  is  a  clear  decrease  in  T;  at  the 
equator  of  approximately  50%  to  a  minimum  value  of 
547K.  This  is  comparable  in  spatial  extent  with  the  BPD 
noted  on  DMSP  and  the  neutral  depletion  on  GRACE. 
Similar  variations  in  T,  occur  during  the  2114  UT  crossing. 

[22]  During  2008  the  plasma  density  in  the  BPD  is 
reduced  to  such  an  extent,  and  the  percentage  of  light  ions 
so  high,  that  accurate  values  of  T,  cannot  be  obtained 
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reliably.  However  up  to  the  point  where  the  uncertainty  in  T; 
becomes  large,  the  temperatures  are  clearly  decreasing.  This 
is  true  of  both  DMSP  and  C/NOFS. 

[23]  Tj  measured  on  DMSP  F16  in  the  BPDs  is  lower 
than  values  reported  during  other  nighttime  solar  minimum 
conditions  when  Te  was  approximately  600-700K  over  the 
Indian  sector  [Bhuyan  et  al.,  2002].  They  are  also  lower 
than  values  predicted  by  the  International  Reference  Iono¬ 
sphere  (IRI)  model  [Gulyaeva  and  Titheridge,  2006]  run  for 
these  periods  (June  2005  and  June  2008)  and  approximate 
location  (330°E  longitude).  These  model  values  are  approx¬ 
imately  1 000K.  We  believe  that  the  discrepancy  between  our 
results  and  past  studies  is  due  to  the  localized  nature  of  BPDs, 
both  in  time  and  space.  Variability  in  ion  and  electron 
densities  and  temperatures  have  been  noted  in  previous 
studies  [Forbes  et  al.,  2000;  Gulyaeva  and  Titheridge, 
2006]  but  no  consistent  study  of  density  variations  as 
functions  of  (1)  solar  cycle,  (2)  season,  and  (3)  location  has 
been  undertaken  until  now. 

5.  Summary 

[24]  We  have  presented  observations  of  BPDs  and  neutral 
density  depletions  which  occur  during  June  solstices,  near 
the  equator  during  solar  minimum  years.  These  observations 
were  made  from  altitudes  of  330  km  to  840  km  using 
several  different  detectors  flown  on  a  number  of  spacecraft. 
The  appearance  of  simultaneous  neutral  density  depletions 
in  the  equatorial  region  narrows  the  range  of  explanations  to 
cooling  or  downwelling  of  the  ionosphere  and  thermosphere 
during  these  times.  This  hypothesis  is  supported  by  direct 
measurements  of  plasma  temperature,  which  also  show 
minima  in  the  equatorial  region,  well  below  the  predicted 
values  based  on  the  IRI  model  or  other  observations.  At 
present  no  mechanism  has  been  proposed  which  accounts 
for  the  basic  formation  of  BPDs  or  their  climatology. 
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[1]  We  report  on  plasma  densities  and  electric  fields 
measured  by  the  C/NOFS  satellite  between  10  and  20  June 
2008.  Midway  through  the  interval,  geomagnetic  conditions 
changed  from  quiescent  to  disturbed  as  a  high  speed  stream 
(HSS)  in  the  solar  wind  passed  Earth.  During  the  HSS 
passage  C/NOFS  encountered  post-midnight  irregularities 
that  ranged  from  strong  equatorial  plasma  bubbles  to 
longitudinally  broad  depletions.  At  the  leading  edge  of  the 
HSS  the  interplanetary  magnetic  field  rapidly  intensified 
and  rotated  causing  auroral  electrojet  currents  to  rise  and 
fall  within  a  few  hours.  As  the  electrojet  relaxed,  C/NOFS 
witnessed  a  rapid  transition  from  a  weakly  to  a  strongly 
disturbed  equatorial  ionosphere  that  lasted  ~10  hours. 
Eastward  polarization  electric  fields  intensified  within  locally 
depleted  flux  tubes.  We  discuss  relative  contributions  of 
gravity-driven  currents,  overshielding  electric  fields  and 
disturbance  dynamos  as  drivers  of  post-midnight  depletions. 
Citation:  Burke,  W.  J.,  O.  de  La  Beaujardiere,  L.  C.  Gentile,  D.  E. 
Hunton,  R.  F.  Pfaff,  P.  A.  Roddy,  Y.-J.  Su,  and  G.  R.  Wilson  (2009), 
C/NOFS  observations  of  plasma  density  and  electric  field 
irregularities  at  post-midnight  local  times,  Geophys.  Res.  Lett., 
36,  L00C09,  doi:10.1029/2009GL038879. 

1.  Introduction 

[2]  Electromechanical  forces  that  couple  the  magnetosphere- 
ionosphere-thermosphere  (MIT)  system  can  cause  intense 
plasma  density  (6N,)  irregularities  to  grow  in  the  low-latitude 
ionosphere.  To  improve  understanding  of  irregularity  forma¬ 
tion,  the  Communication/Navigation  Outage  Forecasting 
System  (C/NOFS)  satellite  was  launched  on  16  April  2008 
into  a  13°  inclined  orbit  with  apogee  and  perigee  at  ^850  and 
400  km,  respectively,  de  La  Beaujardiere  and  C/NOFS 
Science  Definition  Team  [2003]  described  the  mission’s  goals 
and  the  scientific  payload  designed  to  attain  them.  Here  we 
examine  <5N;  variations  sampled  by  the  Planar  Langmuir 
Probe  (PLP)  on  C/NOFS  during  the  solar  minimum  period 
10-20  June  2008  when  perigee  was  on  the  nightside. 

[3]  A  striking,  new  feature  of  C/NOFS  observations  is  the 
detection  of  deep  plasma  depletions  at  topside  altitudes  after 
local  midnight.  Reported  irregularities  appear  as:  (1)  local 
depletions  called  equatorial  plasma  bubbles  (EPBs)  and 

(2)  longitudinally  broad,  depletions  that  we  call  “plasma 
trenches.”  EPBs  are  magnetically  field-aligned,  with  typical 
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east- west  dimensions  of  ~  100  km  |  Basu,  1997]  and  are 
usually  found  at  evening  local  times  (LTs).  Reported  satellite 
detections  of  post-midnight  EPBs  are  relatively  rare  [Burke, 
1979].  Longitudinally  broad  (>500  km)  depletions  develop  in 
the  evening  sector  during  the  main  phases  of  large  magnetic 
storms  when  the  peak  of  the  F-layer  rises  above  the  spacecraft 
[Greenspan  et  al.,  1991],  C/NOFS  has  detected  similar 
phenomena  at  post-midnight  local  times  under  much  less 
stressful  geomagnetic  conditions. 

[4]  EPBs  begin  as  small-amplitude  irregularities  at  bottom- 
side  altitudes  then  intensify  via  a  generalized  Rayleigh-Taylor 
(R-T)  instability  [Balsley  et  al,  1972;  Ott,  1978].  Forces 
driving  the  R-T  instability  include  gravity-driven  currents 
[Eccles,  2004]  and  quasi-DC  electric  fields  (E0)  with  east¬ 
ward  components.  The  main  sources  of  eastward  components 
in  the  dusk  sector  are:  (1)  solar  quiet  (Sq)  current  system 
interactions  with  conductivity  gradients  near  sunset,  and 
(2)  penetration  electric  fields  from  high  latitudes.  Polariza¬ 
tion  electric  fields  (EP)  can  be  much  larger  than  E0,  allow¬ 
ing  bubbles  to  percolate  through  the  topside  at  high  speeds. 
Inside  trenches  E0  may  dominate.  Our  analysis  utilizes 
simultaneous  measurements  of  plasma  densities  by  the 
planar  Langmuir  probe  (PLP)  and  the  Vector  Electric  Field 
Instrument  (VEFI)  [de  La  Beaujardiere  and  C/NOFS  Science 
Definition  Team,  2003]  observed  in  and  near  irregularities. 
VEFI  measures  the  total  electric  field’s  east-west  (Ezonal) 
component  that  includes  E0  and  EP  During  brief  intervals 
when  VEFI’s  pre-amplifiers  were  subject  to  oscillations  data 
outputs  do  not  reflect  geophysical  conditions  and  are  not 
plotted. 

[5]  The  systematics  of  background  Ezonal  at  low  latitudes 
have  been  inferred  from  plasma  flows  detected  by  satellites 
[Fejer  et  al.,  2008]  and  from  incoherent  backscatter  radar 
measurements  [Scherliess  and  Fejer,  1997].  Vertical  drifts 
caused  by  eastward  electric  fields  (Ezonai  >  0)  are  character¬ 
ized  by  pre-reversal  enhancements  (PREs)  near  the  dusk 
terminator.  At  post-sunset  LTs  plasma  initially  rises  in  re¬ 
sponse  to  polarization  charges  [Farley  et  al.,  1986];  and/or 
Sq  current  diversions  [Haerendel  and  Eccles,  1992]  near  the 
terminator.  After  ~20:00  LT  the  direction  of  the  vertical  drift 
usually  reverses,  stabilizing  the  bottomside  F-layer  against 
irregularity  growth.  Ezonal  normally  remains  westward  across 
the  nightside.  Thus  the  PRE  opens  windows  of  opportunity 
for  EPBs  to  form  before  damping  dominates.  PRE  intensities 
have  seasonal,  longitudinal,  and  solar  cycle  dependencies, 
tending  to  be  largest  near  solar-maximum  equinoxes 
[Scherliess  and  Fejer,  1997].  At  solar  minimum  PRE  drifts 
are  weak. 

[6]  Mechanisms  responsible  for  post-midnight  Ezonal  >  0 
at  low  latitudes  include:  (1)  dusk-to-dawn,  overshielding 
electric  fields  that  appear  in  the  recovery  phases  of  storms 
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Figure  1.  Interplanetary  and  auroral  electrojet  activity 
between  1 0  and  20  June  2008,  including:  solar  wind  densities 
(blue)  and  speeds  (red)  measured  near  (top)  Ll5  (middle)  IMF 
BY  (blue)  and  Bz  (red),  and  (bottom)  the  AE  index. 

and  substorms  [Kikuchi  et  al,  2000,  2008]  while  Region  2 
remain  relatively  strong,  and  (2)  disturbance  dynamos  excited 
by  Poynting  and  energetic  particle  fluxes  into  the  auroral 
ionosphere  [Blanc  and  Richmond,  1980;  Fejer  and  Scherliess, 
1995,  1997 ;  Scherliess  and  Fejer,  1997].  Blanc  and  Richmond 
[1980]  argued  that  heating  at  high  latitudes  drives  equatorward 
winds.  Associated  equatorward  Pedersen  currents  accumulate 
polarization  charges  at  mid  latitudes  that  drive  plasma  cir¬ 
culation  in  the  anti-Sq  sense  with  equatorial  electric  fields 
having  eastward  (westward)  components  on  the  nightside 
(dayside).  Fejer  and  Scherliess  [1995]  demonstrated  agree¬ 
ment  between  the  LT  distributions  of  radar  measured  vertical 
drifts  and  the  model  predictions.  C/NOFS  detections  of  post¬ 
midnight  EPBs  provide  opportunities  to  weigh  contributions 
of  these  mechanisms  for  producing  Ezonal  >  0  on  the  night- 
side.  In  principle,  disturbance  dynamo  Ezonal  produces  rela¬ 
tively  sustained  increases  in  post-midnight  EPB  activity 
beginning  a  few  hours  after  the  disturbance  onset  and  lasting 
up  to  30  hours  after  they  abate  [ Scherliess  and  Fejer,  1997], 
Overshielding  fields  occur  early  in  recovery  while  Region  2 
currents  relax  at  slower  rates  than  those  of  Region  1  [Kikuchi 
et  al.,  2000]. 

2.  Observations 

[7]  Figure  1  summarizes  interplanetary  conditions  ob¬ 
served  by  the  Advanced  Composition  Explorer  (ACE)  satel¬ 
lite  and  induced  geomagnetic  responses  during  10-20  June 
2008  (Julian  Days  162-172).  The  top  plot  shows  solar  wind 
densities  Nsw  (blue)  and  speeds  VSw  (red)  from  measure¬ 
ments  by  the  Solar  Wind  Electron  Proton  and  Alpha  Monitor 
[McComas  etal.,  1998].  The  second  shows  the  interplanetary 
magnetic  field  (IMF),  GSM  By  (blue)  and  Bz  (orange) 


components  measured  by  the  Magnetic  Field  Instrument 
[Smith  et  al.,  1998].  The  bottom  plot  is  the  AE  index.  While 
quantities  in  Figure  1  are  hourly  averages,  our  analyses  use 
higher  resolution  data. 

[8]  On  JD  166,  Vsw  rose  from  ^310  km/s  at  16:00  UT  to 
640  km/s  and  remained  >600  km/s  though  JD  170.  Nsw 
data  are  unavailable  prior  to  15:36  UT  on  JD  166.  As  Vgw 
increased,  Nsw  first  rose  to  18.6  then  fell  to  about  2  cm-  . 
Combined  Nsw-Vsw  variations  near  17:00  UT  indicate  that 
a  corotating  interaction  region  (C1R)  formed  at  the  leading 
edge  of  the  HSS  [Balogh  et  al.,  1999;  Tsurutani  et  al., 
2006].  The  IMF  was  weak  both  before  Vsw  increased  and 
inside  the  HSS.  Within  the  CIR  transverse  IMF  components 
reached  ~12  nT  then  rotated  from  northwest  to  south  to 
northwest.  From  JD  162  to  166.5  AE  was  low  44.6  ±  24.3  nT 
with  weak  events  excited  by  brief  southward  excursions.  AE 
responded  strongly  to  IMF  turnings  in  the  CIR  and  HSS.  Peak 
excursions  of  1311  and  1030  nT  occurred  at  21:27  UT  on 
JD  166  and  at  03:55  UT  on  JD  167,  respectively.  It  remained 
high  but  erratic  (237.0  ±  184.5  nT)  as  the  HSS  passed  Earth. 

[9]  C/NOFS  data  are  considered  in  two  stages  with  N;  and 
Ezonal  reported  at  1-s  cadences.  Figure  2  shows  N;  measure¬ 
ments  from  consecutive  orbits  near  the  CIR  impact.  Figure  3 
contains  N;  and  Ezonai  measurements  representative  nightside 
passes  before  and  during  the  HSS  passage.  Unsurprisingly, 
Nj  and  Ezonai  variations  anti-correlate  in  regions  of  plasma 
irregularities.  Before  examining  C/NOFS  data  it  is  useful  to 
consider  orbital  constraints  that  influence  interpretations. 
On  1 0  June  C/NOFS’  perigee  and  ascending  node  were  near 
01:00  and  19:00  LT.  The  ascending  node’s  longitude  pro¬ 
gresses  ~24.2°  to  the  west  on  successive  orbits.  C/NOFS 
thus  provides  3  to  4  snapshots  of  the  nightside  ionosphere  at 
a  given  longitude,  albeit  at  different  latitudes  and  altitudes. 
Orbital  lines  of  apsides  and  nodes  precessed  ~7.5°  to  the  east 
and  ^8.25°  to  the  west  per  day,  respectively.  C/NOFS 
approached  perigee  after  local  midnight.  Locations  of  the 
ascending  node  and  perigee  placed  C/NOFS  north  of  the 
geographic  equator  at  nightside  LTs.  Perigee  was  near 
the  magnetic  equator  in  the  African  to  central  Pacific  sector 
but  at  Appleton  anomaly  magnetic  latitudes  at  eastern 
Pacific  to  South  American  longitudes. 

[10]  Figure  2  shows  N;  measurements  from  late  on 
14  June.  Double  arrows  mark  the  durations  of  C/NOFS 
orbits  878  and  879.  The  letters  M  and  P  mark  midnight 
and  perigee  crossings.  Horizontal  bars  span  periods  when 


UT  20:00  21:00  22:00  23:00  24:00 

Figure  2.  PLP  plasma  densities  measured  between  20:00 
and  24:00  UT  on  14  June  2008.  Times  of  C/NOFS  orbits 
878  and  879  are  indicated  along  with  midnight  meridian  (M) 
and  perigee  (P)  crossings.  Solid  bars  mark  intervals  when 
C/NOFS  was  in  darkness. 
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LT  18:00  20:00  22:00  00:00  02:00  04:00  06:00  08:00 

UT  02:13  02:21  02:29  02:36  02:44  02:52  02:59  03:07 

GLat  -4.1  2.8  8.9  12.4  12.8  9.6  4.0  -2.8 

Gbng  236.7  264.7  292.7  320.8  348.9  17.6  45.0  73.0 

MLat  1.1  12.4  21.4  13.9  -0.2  -3.2  -5.5  -12.6 


LT  18:00  20:00  22:00  00:00  02:00  04:00  06:00  08:00 

UT  03:23  03:32  03:40  03:47  03:55  04:03  04:10  04:17 

GLat  3.4  9.3  12.6  12.6  9.3  3.4  -3.5  -9.4 

Gbng  219.3  247.2  275.2  303.3  331.4  359.5  27.5  55.5 

MLat  6.5  15.8  23.9  22.6  4.6  -9.6  -14.9  -19.3 


Figure  3.  (top)  Zonal  electric  fields  and  (middle)  plasma  densities  measured  by  C/NOFS  during  orbits  (a)  822  and  (b)  882, 
plotted  as  functions  of  LT,  UT,  geographic  latitude,  geographic  longitude,  and  geomagnetic  latitude.  Orbit  822  illustrates 
predisturbance  conditions;  orbit  882  occurred  shortly  after  the  C1R  passage,  (bottom)  Traces  indicate  geographic  altitudes 
(solid  lines)  of  C/NOFS  and  the  apex  altitudes  of  magnetic  field  lines  on  which  measurements  were  made. 


C/NOFS  was  in  darkness.  During  the  eclipsed  portions  of 
orbits  878  and  879  C/NOFS  crossed  longitudes  —10°  to 
105°E  and  —35°  to  80°E,  respectively.  Thus,  both  passes 
covered  African  to  Indian  Ocean  longitudes.  Except  for 
one  depletion  near  21:20  UT  (00:45  LT,  ~50°E)  bN, 
amplitudes  were  small  during  orbit  878.  Conversely,  on 
orbit  879  C/NOFS  encountered  continuous  large  /NpN;  > 
100  irregularities  from  eastern  Atlantic  through  African 
longitudes.  Pre-midnight  EPBs  apparent  on  orbit  879  are 
singular  exceptions  to  the  post-midnight  rule  for  these  ten 
days.  Expanded  views  of  Nt  (not  given)  indicate  that  local 
bNi  depressions  were  110  to  480  km  wide.  The  African 
sector  is  common  to  both  observations.  Thus,  differences 
observed  between  orbit  878  and  879  are  time  dependent. 
Note  that  the  first  AE  peak  occurred  ~1  hour  before  the 
start  of  orbit  879. 

[1 1]  Figure  3  shows  Ezonal  (Figure  3,  top)  and  N; 
(Figure  3,  middle)  sampled  at  similar  longitudes  during 
C/NOFS  orbits  822  (Figure  3a)  and  882  (Figure  3b).  These 
data  were  acquired  during  the  magnetically  quiet/disturbed 
JD  163/167,  respectively  and  are  plotted  as  functions  of 
LT  beginning  at  ascending  nodes,  UT,  geographic  latitude 
(GLat),  geographic  longitude  (GLong)  and  corrected  mag¬ 
netic  latitude  (MLat).  The  bottom  plots  indicate  spacecraft 
altitudes  (solid  lines)  and  apex  altitudes  of  magnetic  field 
lines  (dashed  lines).  Traces  converge  where  C/NOFS  was 
near  the  magnetic  equator. 

[12]  Plasma  densities  shown  in  Figure  3a  increased  at 
pre-midnight  LTs  as  the  C/NOFS  altitude  decreased.  Even¬ 
ing  sector  irregularities  appear  as  weak  local  increases  and 
decreases.  Starting  at  ~23:00  LT  while  C/NOFS  was 


descending,  N,  decreased  gradually.  After  local  midnight 
N,  fell  precipitously  near  perigee.  From  ~01:00  LT  through 
dawn  C/NOFS  crossed  a  series  of  (factors  of  10  to  100) 
EPB  depletions;  Ezonal  was  irregular  and  eastward  within  the 
depletions.  C/NOFS  remained  close  to  the  magnetic  equator 
as  its  trajectory  rose  from  perigee.  Figure  3b  (middle)  shows 
a  pre-midnight  N;  enhancement  followed  by  post-midnight 
irregularities  and  depletions,  but  EPB  signatures  were 
infrequent.  Rather,  C/NOFS  crossed  a  wide  (~3500  km) 
plasma  density  trench  from  ^03:53  to  04:01  UT  (01:45- 
04:00  LT)  on  JD  167.  Across  this  deep  depletion,  the 
average  Ezonal  was  ~2.4  mV/m  eastward  (Figure  3b,  top). 
After  03:56:15  UT,  (~02:30  LT)  when  the  plasma  density 
abruptly  fell  to  ~20  cm  3,  Ezonal  variations  increased  dra¬ 
matically  but  with  no  correspondence  to  simultaneous 
changes  in  measured  <5Ni.  At  these  low  densities  the  Debye 
length  grew  to  ~1  m.  Thus,  observed  Ezonal  fluctuations 
probably  reflect  responses  to  differences  in  VEFl’s  sheath 
impedances.  We  note,  but  do  not  show,  that  the  large- 
amplitude  irregularities  of  orbit  879  appear  to  have  evolved 
into  the  broad  trench.  During  nightside  parts  of  orbits  880 
and  881,  C/NOFS  encounters  with  topside  plasma  dimin¬ 
ished  and  the  longitudinal  scale  sizes  of  N;  depletions 
widened.  The  deep  depletion  of  orbit  882  narrowed  during 
orbit  883  then  disappeared  as  the  affected  ionosphere  rotated 
into  the  dayside. 

[13]  Data  acquired  during  nightside  passes  of  Figure  3 
also  illustrate  relationships  between  N;  (blue)  and  Ezonal  (red). 
During  orbit  822  bN;  and  <5Ezonal  strongly  anti-correlated, 
Ezonal  intensified  within  N;  depletions.  This  is  consistent  with 
C/NOFS’  passing  through  upwelling  EPBs  permeated  by 
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strong  polarization  electric  fields  [Ott,  1978].  This  relation¬ 
ship  continued  during  orbit  882,  from  03:54-03:55  UT, 
while  C/NOFS  was  in  the  western  part  of  the  deep  plasma 
trench.  However,  at  ^03:56:13  UT  as  N,  rapidly  decreased 
from  to  ~15  cm-3,  />Ezonal  underwent  large  amplitude 
oscillations  related  to  large  sheaths  around  VEF1.  The  quasi 
DC  portion  of  Ezonal  was  ~2.4  mV/m  across  the  orbit  882 
plasma  trench,  consistent  with  the  local  F  layer  being  lifted 
above  C/NOFS. 


3.  Discussion 

[14]  This  study  focused  on  dynamics  of  the  nightside, 
equatorial  ionosphere  before  and  during  passage  of  a  HSS 
in  the  solar  wind.  During  this  period  C/NOFS’  perigee  was 
at  post-midnight  LTs.  Interplanetary  and  AE  data  show  three 
distinct  intervals,  (1)  geomagnetic  quiet,  JD  162-166.5, 
(2)  responses  to  the  C1R  magnetic  field  rotation  JD  166.5- 
167,  and  (3)  HSS  induced  activity  JD  167-170. 

[15]  Figure  3a  (middle)  shows  that  from  the  western 
Atlantic  across  Africa  to  the  central  Pacific  (~  1 60°E)  C/NOFS 
detected  pre-disturbance  EPBs  and  large  depletions  after 
midnight.  At  other  longitudes  C/NOFS  was  away  from  the 
magnetic  equator.  Even  during  sustained  periods  of  magnetic 
quiet  in  solar  minimum  bottomside  irregularities  grew  into 
EPBs.  This  raises  a  question.  What  mechanism  supports  the 
development  of  EPBs  at  post-midnight  LTs  during  quiet 
intervals  when  neither  overshielding  nor  disturbance  dynamo 
Ezonal  fields  can  operate?  Solar  UV  fluxes  driving  dayside 
dynamos  and  consequent  PRE  signatures  were  also  weak. 
With  no  post-dusk  westward  Ezonal  to  stabilize  the  bottomside 
ionosphere,  the  gravitational  drift  current jG  =  njrqg  x  B)/B2 
has  no  competitors.  Normally  insufficient  tune  is  available  for 
jo  to  drive  bottomside  irregularities  into  the  nonlinear  regime 
before  Ezonal  turns  westward  [. Basu ,  1997].  If  Ezonal «  0,  data 
suggest  that  enough  time  becomes  available  for  EPBs  to 
grow  and  reach  altitudes  >400  km  after  midnight. 

[16]  Figure  2  shows  that  the  nightside  ionosphere 
responded  to  the  CIR  passage.  Except  for  the  high  number 
of  EPBs  encountered  during  orbit  879  the  upper  envelopes  of 
Ni  traces  for  orbits  878  and  879  were  similar.  The  R-T  growth 
rate  equation  [Ott,  1978]  suggests  that  an  eastward  Ezonal  was 
present  on  the  nightside  during  the  hour  before  C/NOFS  orbit 
879.  AE  first  peaked  at  21:27  UT  then  decreased  at  a  rate  of 
283  nT/hr.  Scherliess  and Fejer  [1997]  estimated  that  several 
hours  are  needed  for  disturbance-dynamo  effects  to  reach 
the  equator.  On  the  other  hand,  overshielding  electric  fields 
appear  soon  after  the  IMF  turns  northward  [Kikuchi  et  al., 
2000].  We  suggest  that  overshielding  was  primarily  respon¬ 
sible  for  changes  seen  between  orbits  878  and  879. 

[17]  During  orbit  882  C/NOFS  crossed  a  broad  and  deep 
depletion  with  average  Ezonal  rss  2.4  mV/m.  With  an  equato¬ 
rial  magnetic  field  of  ~3 1 ,000  nT,  this  corresponds  to  upward 
drifts  of  ^77  m/s.  It  would  take  ~22  minutes  for  the  F  layer  to 
rise  100  km.  Thus,  within  this  deep  trench  C/NOFS  sampled 
plasma  that  originated  far  below  the  peak  height  of  the 
quiescent  F  layer.  Between  JD  167  and  170  C/NOFS  wit¬ 
nessed  several  cycles  of  plasma  and  AE  quieting  and  activa¬ 
tion.  With  available  information  it  is  difficult  to  distinguish 
between  overshielding  and  disturbance-dynamo  contribu¬ 
tions.  However,  the  average  Ezonal  «  2.4  mV/m  measured 
across  the  deep  plasma  trench  seems  to  be  a  disturbance 


dynamo  effect.  Early  on  JD  167  the  downward  trend  of  AE 
reversed,  reaching  a  second  peak  of  1030  nT  at  03:55  UT.  It 
appears  unlikely  that  an  overshielding  Ezonai  operated  in 
the  hour  prior  to  orbit  882.  Within  the  bounds  of  present 
understanding  only  the  disturbance  dynamo  can  explain  the 
strong  eastward  Ezonal  measured  during  orbit  882. 

4.  Conclusions 

[is]  This  study  of  C/NOFS  measurements  indicates  that 
during  June  2008  EPBs  and  deep  plasma  trenches  formed  at 
post-midnight  LTs.  These  features  were  present  during  an 
extended  period  of  magnetic  quiet,  but  intensified  after  the 
CIR’s  coupling  to  the  magnetosphere-ionosphere  weakened 
and  allowed  the  auroral  electrojet  to  relax.  Observed  EPB 
activity  waxed  and  waned  as  the  HSS  passed  Earth.  Our 
analysis  supports  four  conclusions: 

[19]  1.  Active  post-midnight  EPBs  carry  eastward  polar¬ 
ization  electric  fields  and  indicate  that  irregularities  can 
form  and  reach  C/NOFS  altitudes  even  during  magnetically 
quiet  intervals  in  solar  minimum.  Under  these  conditions 
nim,(g  x  B)  currents  acquire  sufficient  time  to  drive 
nonlinear  R-T  instabilities. 

[20]  2.  Intensification  of  EPB  activity,  observed  about  an 
hour  after  the  CIR-induced  peak  in  AE,  was  probably  caused 
by  eastward  overshielding  electric  fields. 

[21]  3.  The  plasma  trench  of  orbit  882  with  Ezonai  « 
2.4  mV/m  was  observed  while  AE  was  increasing  and 
thus  appears  to  be  mostly  a  disturbance  dynamo  effect. 

[22]  4.  Continued,  albeit  sporadic,  detections  of  EPBs  and 
plasma  trenches  during  the  four  days  of  the  HSS  passage 
suggest  that  both  overshielding  and  disturbance  dynamos 
remained  active.  Presently  available  information  is  insuffi¬ 
cient  to  distinguish  their  distinctive  contributions. 
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[l]  Using  electric  field  measurements  as  inputs,  the 
assimilative  physics-based  ionospheric  model  (PBMOD) 
successfully  reproduced  density  depletions  observed  at  early 
morning  local  times  during  four  consecutive  orbits  of  the 
Communication/Navigation  Outage  Forecasting  System 
(C/NOFS)  satellite  on  17  June  2008.  However,  the  PBMOD 
running  with  plasma  drift  data  from  empirical  models  as 
inputs  predicted  neither  plasma  depletions  nor  irregularities 
on  this  day.  Coincident  over  flights  of  a  large  depletion  by 
C/NOFS  and  the  DMSP-F17  satellite  allow  estimates  of  its 
longitudinal  and  latitudinal  scale  sizes.  The  satellite-based 
estimates  are  shown  to  be  in  reasonable  agreement  with 
PBMOD  predictions.  The  model’s  reproduction  of  observed 
temporal  and  spatial  distributions  of  plasma  depletions 
suggests  that  our  assimilative  technique  can  be  used  to 
enhance  space-weather  forecasts.  Citation:  Su,  Y.-J.,  J.  M. 
Retterer,  O.  de  La  Beaujardiere,  W.  J.  Burke,  P.  A.  Roddy,  R.  F. 
Pfaff  Jr.,  G.  R.  Wilson,  and  D.  E.  Hunton  (2009),  Assimilative 
modeling  of  equatorial  plasma  depletions  observed  by  C/NOFS, 
Geophys.  Res.  Lett.,  36,  L00C02,  doi:10.1029/2009GL038946. 

1.  Introduction 

[2]  Properties  of  low-latitude  plasma  irregularities,  com¬ 
monly  observed  at  post-sunset  local  times,  have  been 
studied  extensively  using  ground  and  satellite-based  mea¬ 
surements  [e.g.,  Fejer  and  Kelley,  1980;  Fejer,  1997].  The 
pre-reversal  enhancement  of  upward  plasma  drifts,  due  to 
eastward  polarization  electric  fields,  is  the  primary  cause  of 
these  irregularities  [e.g.,  Richmond,  1995].  On  the  other 
hand,  reported  plasma  depletions  at  post-midnight  and  early 
morning  local  times  are  scarce.  The  majority  of  reported 
events  occurred  during  storm  time  conditions  [e.g.,  Burke, 
1979;  Su  et  al.,  2004].  In  a  rare  occurrence  during  solar 
maximum,  ROCSAT-1  simultaneously  observed  an  intense 
reduction  in  ion  density  with  super-cool  temperatures, 
accompanied  by  large  downward,  eastward,  and  field- 
aligned  flows  near  the  dawn  meridian  [Su  et  al.,  2004]. 
Burke  [1979]  suggested  that  dawn  bubbles  may  be  caused 
by  substorm-generated  eastward  electric  fields  at  the  equator 
or  by  travelling  ionospheric  disturbances.  In  a  theoretical/ 
numerical  study,  Eccles  [2004]  showed  a  persistent  10  m  s_1 
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pre-sunrise  uplift  of  the  ionosphere  due  to  gravity-driven 
current.  In  a  statistical  study  by  Sobral  et  al.  [1999],  two 
cases  were  found  between  1980  and  1992  in  the  post¬ 
midnight  sector  during  quiet  times  in  which  plasma  deple¬ 
tions  were  observed.  C/NOFS  frequently  encounters  deep 
plasma  depletions  at  the  dawn  terminator,  particularly  in 
June  and  September  2008.  We  have  analyzed  density  and 
electric  field  irregularities  at  post-midnight  local  times 
during  a  10-day  period  in  June  2008  during  which  a  fast 
stream  in  the  solar  wind  passed  the  Advanced  Composition 
Explorer  (ACE)  satellite.  This  paper  describes  equatorial 
ionospheric  responses  observed  in  the  early  morning  local 
times  on  17  June  2008  while  a  solar  wind  fast  stream  was 
passing  Earth  for  comparison  with  predictions  of  the  physics- 
based  ionospheric  model  (PBMOD)  [Retterer,  2005;  Retterer 
et  al.,  2005].  PBMOD  operates  in  both  assimilative  and 
climatological  modes.  The  former  mode  is  done  with  electric 
field  measurements  from  the  Vector  Electric  Field  Instrument 
(VEFI)  on  C/NOFS.  The  latter  has  no  electric-field  informa¬ 
tion  beyond  that  provided  by  climatological  models  [e.g., 
Scherliess  and  Fejer,  1999].  Plasma  densities  obtained  from 
both  types  of  simulation  are  compared  with  measurements 
by  the  Planar  Langmuir  Probe  (PLP)  on  C/NOFS. 

2.  Model  Description 

[3]  The  simulation  tool  used  in  this  study  was  developed 
by  Retterer  [2005]  over  the  last  decade  at  the  Air  Force 
Research  Laboratory  to  support  the  C/NOFS  mission  anal¬ 
ysis.  This  4-D  model  of  the  low-latitude  ionosphere  solves 
the  continuity  equation  as  functions  of  position  and  time 
along  magnetic  field  lines  over  a  range  of  apex  altitudes. 
The  climatological  version  of  the  simulation  is  carried  out 
with  initial  conditions  such  as  neutral  winds  and  temper¬ 
atures,  ion  production  and  loss  rates,  and  plasma  density 
and  temperature  prescribed  by  various  empirical  models 
[Retterer  et  al.,  2005].  The  assimilative  version  input  data 
consist  of  available  satellite  and/or  ground  measurements  as 
inputs.  In  this  study,  the  only  ingested  data  are  E  x  B  drift 
velocities  calculated  from  VEFI  electric  field  measure¬ 
ments  and  the  International  Geomagnetic  Reference  Field 
(IGRF)  model  [Mans  et  ah,  2005],  We  note  that  the  VEFI 
instrument  also  provides  magnetic  field  measurements.  The 
magnetic  field  strength  obtained  from  IGRF  is  within  1% 
of  the  measured  value.  Whenever  electric  field  data  were 
unavailable,  such  as  simulation  regions  far  away  from  the 
spacecraft  trajectory,  we  used  climatological  values  from 
an  empirical  drift  model  based  on  Jicamarca  radar  and  the 
Atmospheric  Explorer  satellite  data  [Scherliess  and  Fejer, 
1999].  Retterer  [2005]  provides  a  fuller  description  of  the 
assimilative  procedures. 
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Figure  1.  Four  one-hour  segments  of  data  at  (a)  04:00-05:00,  (b)  05:42-06:42,  (c)  07:30-08:30,  and  (d)  09:00-10:00  UT. 
From  top  to  bottom  plots:  ion  densities,  vertical  and  zonal  drift  velocities,  and  satellite  altitude,  latitude,  and  longitude. 


[4]  Although  PBMOD  was  designed  as  a  tool  to  study 
both  coarse-scale  global  ionospheric  structures  and  fine- 
scale  localized  plasma  turbulences,  this  paper  focuses  only 
on  large-scale  density  depletions  observed  at  post-midnight 
local  times  on  17  June  2008  (day  of  year  169).  While  the 
average  Kp  index  was  ~2  on  this  solar-minimum  day,  the 
auroral  electroject  (AE)  index  showed  large  amplitude 
variations  typical  of  fast  stream  intervals  [ Tsurutani  and 
Gonzalez,  1987].  The  simulated  low-latitude  ionosphere  is 
divided  into  36  longitudinal  sectors  of  10°  width.  The 
latitudinal  grid  has  201  steps  along  each  field  line.  The 
vertical  grid  has  60  points  varying  in  size  and  ranges  for 
field-line  apex  altitudes  between  140  and  4000  km.  The 
simulation  time  step  is  6  min.  Throughout  this  report  sim¬ 
ulation  results  and  satellite  observations  are  presented  in 
geographic  coordinates. 

3.  Results 

[5]  Although  PBMOD  was  run  in  both  the  assimilative 
and  climatological  modes  for  the  entire  day  of  17  June 
2008,  here  we  focus  on  four  one-hour  segments  that  include 
consecutive  C/NOFS  nightside  passes  labeled  912-915  in 
Figures  la- Id.  Black,  blue,  and  red  lines  in  Figures  la- Id 
represent  C/NOFS  observations,  climatology  results,  and 
assimilative  model  outputs,  respectively.  Ion  densities  are 
plotted  followed  by  the  vertical  and  zonal  drift  velocities; 
positive  values  indicate  upward  and  eastward  drifts.  Before 
assimilating  electric  fields  into  the  global  model,  we  applied 
a  low-pass  Fourier  filter  [Retterer,  2005]  with  a  5-min 
window  to  1 -second  averaged  VEF1  and  PLP  data  to 
remove  small-scale  variations  (represented  by  the  black 
lines  in  top  three  plots  of  Figure  1).  The  vertical  dashed 
lines  in  the  top  three  plots  of  Figure  1  mark  times  when 
C/NOFS  crossed  the  06:00  LT  meridian.  The  C/NOFS 
altitude  and  the  apex  altitude  of  the  field  line  it  occupies 
are  represented  by  black  and  red  lines  in  the  fourth  plot 
from  the  top.  The  satellite  latitude  and  the  magnetic  equator 
are  plotted  as  black  and  red  lines  in  the  fifth  panel. 


Intersections  of  the  red  and  black  lines  in  plots  4  and  5  of 
Figure  1  mark  satellite  crossings  of  the  magnetic  equator. 
Black  and  red  lines  in  Figure  1  (bottom)  indicate  the 
longitudes  and  local  times  of  the  spacecraft;  the  horizontal 
dotted  lines  denote  06:00  and  18:00  LT. 

[6]  As  seen  in  Figure  1  (top)  large-scale  density  deple¬ 
tions  were  observed  during  four  consecutive  satellite  orbits 
when  eastward  electric  fields  (e.g.,  upward  velocities)  were 
observed  by  VEFI  (black  lines  in  the  second  plot  of 
Figure  1)  when  or  a  few  minutes  before  the  satellite  encoun¬ 
tered  the  depletion.  C/NOFS  crossed  all  of  the  depletions 
presented  here  at  early  morning  local  times  while  flying  at 
or  near  the  magnetic  equator.  PBMOD  reproduced  the  main 
qualitative  features  of  the  pre-dawn  sector  depletions  (red 
lines)  observed  by  PLP  (black  lines)  based  on  assimilations 
of  available  electric  field  data  obtained  from  VEFI  (red  lines 
in  the  2nd  and  3rd  plots  of  Figure  1).  Conversely,  no  plasma 
depletions  or  irregularities  were  predicted  while  exercising 
the  climatological  version  (blue  lines)  with  input  drift  veloc¬ 
ities  derived  from  empirical  models. 

[7]  It  is  impossible  to  specify  both  the  latitude  and 
longitude  scale  sizes  of  plasma  depletions  with  measure¬ 
ments  from  a  single  satellite.  Fortunately,  the  FI 7  satellite  of 
the  Defense  Meteorological  Satellite  Program  (DMSP) 
crossed  the  same  depletion  structure  at  a  higher  altitude 
(860  km)  ~  10  min  prior  to  the  C/NOFS  encounter  at 
420  km  (cf.  Figure  Id).  The  tracks  of  the  two  satellite  orbits 
are  illustrated  in  Figure  2b.  DMSP-F17  crossed  the  depletion 
moving  north  to  south  at  09:12-09:36  UT,  while  C/NOFS 
traversed  it  moving  west  to  east  at  09:24-09:48  UT. 
Corresponding  density  measurements  at  1  sec  resolution 
from  DMSP  and  C/NOFS  are  plotted  in  Figures  2a  and  2c, 
respectively.  Note  that  the  time  scale  increases  downward  in 
Figure  2a.  The  dotted  lines  in  Figures  2a  and  2c  mark  our 
estimates  of  the  depletion  boundaries.  If  we  assume  that  the 
plasma  density  structure  did  not  change  significantly  during 
the  10  min  time  offset,  the  dimensions  of  the  depletion  are 
about  14°  in  longitude  and  21°  in  latitude,  corresponding  to 
the  heavy-line  segments  of  the  trajectories  in  Figure  2b. 
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[8]  A  comparison  of  simulation  results  with  the  two 
satellite  observations  is  presented  in  Figure  3.  It  shows  0+ 
density  outputs  from  PBMOD  calculated  using  assimilated 
VEF1  (Figure  3a)  and  climatological  values  for  the  velocity 
drifts  (Figure  3b).  Recall  that  each  time  step  of  PBMOD 
simulations  is  6  min.  As  indicated  above  Figure  3,  the  2D 
plots  represent  conditions  at  09:36  UT  and  4°S.  Red  stars 
indicate  C/NOFS’  location  at  9:35  UT  when  it  sampled  the 
density  minimum  (295°E,  5°S);  black  stars  mark  DMSP- 
F17’s  location  at  09:25UT  (300°E,  5°S).  Climatology-based 
model  results,  shown  in  Figure  3b,  indicate  low  nightside 
densities  in  the  bottom-side  ionosphere.  Just  prior  to  sunrise 
the  predicted  altitude  of  the  peak  of  the  F2  layer  ( hmF2 )  is  at 
~300  km.  By  contrast,  assimilation-based  results  displayed 
in  Figure  3a  indicate  that  the  ionosphere  was  lifted  upward. 
Near  times  when  the  two  satellites  encountered  the  deep 
depletion  hmF2  had  reached  ~900  km  altitude.  DMSP-F17 
was  flying  below  the  F2  peak.  The  longitudinal  extent  of  the 
depletion  is  ~30°.  In  order  to  demonstrate  the  predicted 
latitudinal  dimension  of  the  depletion,  Figure  3  c  shows 
assimilation  (red)  and  climatological  (black)  predictions  of 
peak  F-layer  densities  ( NmF2 )  plotted  as  functions  of  latitude 
at  09:36  UT  and  285°E.  The  climatologically  predicted 
latitudinal  separation  between  the  two  NmF2  maxima  (the 
equatorial  ionospheric  anomaly)  was  18°.  The  latitudinal 
width  of  density  valley  increased  to  32°  in  assimilation- 
model  calculations.  The  size  of  the  density  depletion  esti¬ 
mated  from  the  two  satellites  shown  in  Figure  2  is  within  the 
dimension  of  the  depletion  structure  derived  from  the 
assimilative  simulation  results  shown  in  Figures  3a  and  3c. 
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Figure  2.  (a)  Total  ion  densities  measured  by  SSIES  on 
board  DMSP-F17  from  9:12  to  9:36  UT,  where  the  y-axis 
label  include  UT,  longitude,  latitude,  and  altitude  information 
of  the  satellite,  (b)  DMSP-F17  crossed  the  equator  from 
north  to  south,  while  C/NOFS  traversed  South  America 
from  west  to  east  10  min  after  the  DMSP  pass,  (c)  The  ion 
densities  observed  by  PLP  on  board  C/NOFS  from  9:24  to 
9:48  UT.  The  boldfaced  cross  in  Figure  2b  indicates  an 
estimated  depletion  region  corresponding  to  the  dotted  lines 
in  Figures  2a  and  2c. 
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Figure  3.  (a)  2D  density  map  obtained  from  the  assimila¬ 
tion  model,  where  the  x-  and  y-  axis  represent  longitude  and 
altitude,  respectively,  (b)  2D  density  map  based  on  the 
climatology  model  output.  The  black  and  red  star  symbols 
in  Figures  3a  and  3b  denote  the  DMSP-F17  and  C/NOFS 
satellite  locations  when  the  minimum  densities  were 
observed  at  9:25  and  9:35  UT,  respectively.  The  black  dots 
indicate  local  hmF2.  (c)  ID  plot  of  NmF2  as  a  function  of 
geographical  latitude,  where  the  red  and  black  lines  represent 
results  from  the  assimilation  and  climatology  models, 
respectively. 


[9]  An  examination  of  the  temporal  evolution  of  density 
maps  derived  from  PBMOD  (see  Animation  SI  of  the 
auxiliary  material)  indicates  that  the  density  depletions 
shown  in  Figures  la- Id  should  be  regarded  as  passing 
through  the  same  structure  during  four  consecutive  C/NOFS 
orbits. 1  This  structure  expands  or  shrinks  in  response  to  the 
history  of  vertical  plasma  drifts  (zonal  electric  fields). 
Figure  3  a  shows  that,  at  longitudes  >  310°E,  the  bottom- 
side  of  the  F-layer  began  to  refill  with  plasmas  created  after 
sunrise,  while  topside  densities  remained  low.  Several 
minutes  later,  the  bottom-side  density  depletion  (<400  km 
between  280°E  and  310°E  longitude)  extended  further  into 
normally  topside  altitudes  due  to  large  upward  velocities 
observed  at  09:36  UT  (second  plot  of  Figure  Id).  However, 


'Auxiliary  materials  are  available  in  the  HTML,  doi:  10.1 029/ 
2009GL038946. 
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the  longitudinal  width  of  the  depletion  continued  to  decrease 
with  time  as  its  most  eastward  portion  crossed  the  dawn 
terminator  where  the  ionosphere  was  replenished  with  new 
plasmas  generated  by  photoionization. 

4.  Summary  and  Discussion 

[10]  The  assimilative  version  of  PBMOD  that  had  VEFI 
inputs  reproduced  qualitative  large-scale  features  of  density 
depletions  observed  near  dawn  by  the  PLP  during  four 
consecutive  C/NOFS  orbits  on  17  June  2008.  No  plasma 
depletions  or  irregularities  were  predicted  by  the  climatology 
version  of  the  model.  DMSP-F17  and  C/NOFS  made  north- 
to-south  and  west-to-east  cuts  through  the  same  density 
depletion  structure  within  1 0  min  of  each  other.  At  the  times 
of  these  encounters  the  latitudinal  and  longitudinal  dimen¬ 
sions  of  the  depletion  are  ~21°  and  14°,  respectively.  This 
observational  result  is  within  the  estimated  32°  by  30°  size 
of  the  density  depletion  derived  from  assimilative  modeling 
results.  Consistent  with  observations  of  the  depletion  at 
420  and  860  km,  the  assimilation  indicated  that  hmF2 
reached  altitudes  near  900  km. 

[n]  A  significant  advantage  of  the  4D  physics-based 
modeling  is  to  provide  reasonable  global  information  of 
ionospheric  responses  when  and  where  in-situ  satellite 
observations  are  not  available.  Assimilative  modeling  results 
confirm  that  C/NOFS  sampled  the  same  density  depletion  in 
four  orbits  between  04:00  and  10:00  UT.  As  the  depletion 
crossed  the  dawn  terminator  it  began  to  refill  with  fresh 
plasmas  from  the  lower  ionosphere.  The  depletion  structure 
expanded  and  contracted  in  response  to  the  history  of  zonal 
electric  fields.  The  uplift  of  the  ionosphere  is  due  to  an 
eastward  electric  field,  while  westward  electric  fields  act  to 
stabilize  the  nightside  ionosphere.  The  cause  of  the  large- 
scale  dawn  density  depletion  is  similar  to  that  of  post-sunset 
plasma  irregularities  or  bubbles,  however,  their  generation 
mechanisms  are  different. 

[12]  Plasma  density  depletions  discussed  in  this  paper 
were  observed  during  solar  minimum  when  the  average 
F10.7  was  ~67.  The  gravity-driven  current  suggested  by 
Eccles  [2004]  may  be  the  dominant  effect  to  explain  the 
dawn  depletion  during  extended  quiet  geomagnetic  condi¬ 
tions.  Upward  drifts  of  10-15  m  s-1  at  dawn  due  to 
niig  x  B  currents,  as  estimated  by  Eccles  [2004],  are 
considerably  smaller  than  those  inferred  from  VEFI  mea¬ 
surements  on  17  June  2008.  We  cannot  rule  out  contributions 
from  disturbance  dynamos  [ Scherliess  and  Fejer,  1997] 
excited  after  a  corotating  interaction  region  (C1R)  at  the 
leading  edge  of  a  fast  solar  wind  stream  observed  by  ACE  at 
~16  UT  on  14  June.  Large  excursions  of  the  AE  index 
occurred  at  2 1  UT  on  1 6  June  and  at  2  UT  on  1 7  June.  Over¬ 
shielding  eastward  fields  may  be  induced  when  Region  2 
currents  exceed  those  of  Region  1  when  AE  decreases 
rapidly  [ Kikuchi  et  al.,  2000].  Further  investigation  is 
needed  in  order  to  separate  various  contributions  in  the 
generation  of  eastward  electric  fields. 

[13]  Small-scale  variations  in  the  drift  velocities  have  been 
smoothed  by  a  low-pass  filter  prior  to  use  as  input  for  the 
global  run  of  PBMOD  in  order  to  maintain  a  reasonable 
model  run-time  for  forecasting  purpose.  As  a  result,  this 
study  is  limited  in  scope  to  an  examination  of  large  density 
structures  in  the  ambient  ionosphere.  The  spatial  and  tem¬ 


poral  resolution  of  the  global  PBMOD  run  can  be  enhanced 
with  increased  computing  power.  The  small-scale  density 
irregularities,  such  as  equatorial  bubbles  or  spread  F,  can  be 
further  investigated  by  the  regional  plume  model  algorithms 
in  PBMOD  [Retterer,  2005], 

[14]  The  neutral  wind  plays  an  important  role  in  the 
thermosphere-ionosphere  dynamics.  In  future  studies,  we 
plan  to  incorporate  measurements  obtained  from  the  Neutral 
Wind  Meter  on  C/NOFS  for  PBMOD  assimilation  to  inves¬ 
tigate  the  neutral  wind  effect  and  to  further  validate  the  space 
weather  forecasting  capabilities. 
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2301SDA5,  Air  Force  contract  FA8718-08-C-0012  with  Boston  College, 
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[i]  The  Communication/Navigation  Outage  Forecasting 
System  (C/NOFS)  satellite  was  launched  in  2008,  during 
solar  minimum  conditions.  An  unexpected  feature  in  the 
C/NOFS  plasma  density  data  is  the  presence  of  deep  plasma 
depletions  observed  at  sunrise  at  all  satellite  altitudes. 
Ionospheric  irregularities  are  often  embedded  within  these 
dawn  depletions.  Their  frequencies  strongly  depend  on 
longitude  and  season.  Dawn  depletions  are  also  observed  in 
coincident  satellite  passes  such  as  DMSP  and  CHAMP.  In 
one  example  the  depletion  extended  50°  x  14°  in  the  N-S 
and  E-W  directions,  respectively.  These  depletions  are 
caused  by  upward  plasma  drifts  observed  in  C/NOFS  and 
ground-based  measurements.  The  reason  for  these  upward 
drifts  is  still  unresolved.  We  discuss  the  roles  of  dynamo 
electric  fields,  over-shielding,  and  tidal  effects  as  sources 
for  the  reported  depletions.  Citation:  de  La  Beaujardiere,  O., 
et  al.  (2009),  C/NOFS  observations  of  deep  plasma  depletions  at 
dawn,  Geophys.  Res.  Lett.,  36,  L00C06,  doi:10.1029/ 
2009GL038884. 

1.  Introduction 

[2]  The  C/NOFS  satellite  was  launched  in  April  2008  as 
an  Air  Force  mission  to  forecast  ambient  plasma  densities 
and  irregularities  in  the  equatorial  ionosphere  \de  La 
Beaujardiere  et  al.,  2004].  Instruments  on  the  satellite 
measure  electric  fields,  plasma  characteristics,  neutral 
winds,  and  the  strength  of  scintillation-producing  irregular¬ 
ities.  The  early  phase  of  C/NOFS  observations  corresponds 
to  the  solar  cycle  minimum  (FI 0.7  =  68).  The  pre-reversal 
enhancement  in  the  vertical  plasma  drift,  responsible  for 
early-evening  irregularities,  rarely  appears  at  such  times. 
Instead,  there  is  a  high  occurrence-rate  of  irregularities  after 
midnight.  In  addition,  distinctive  deep  plasma  depletions  are 
also  observed  at  sunrise.  This  paper  describes  and  analyzes 
these  unexpected  dawn  depletions. 

2.  Observations  of  Dawn  Depletions 

[3]  This  paper  focuses  on  data  acquired  by  two  C/NOFS 
sensors.  The  Planar  Langmuir  Probe  (PLP)  measures  the 
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ambient  ionospheric  densities  and  electron  temperatures.  It 
also  provides  high  temporal  (512  Hz)  and  spatial  (~13  m) 
resolution  measurements  of  density  irregularities.  The  Vec¬ 
tor  Electric  Field  Instrument  (VEFI)  measures  the  AC  and 
DC  electric  and  magnetic  fields.  Electric  fields  are  measured 
with  three  orthogonal  20  m  tip-to-tip  booms.  Initial  compar¬ 
isons  between  VEFI  and  PLP  with  observations  from  the 
Jicamarca  incoherent  scatter  radar  (ISR)  showed  good 
agreement,  supporting  their  validity. 

[4]  Figure  1  provides  our  first  example  of  a  deep  plasma 
depletion  observed  at  sunrise  during  C/NOFS  orbit  #915,  on 
June  8,  2008  at  09:34  UT,  while  the  satellite  was  at  410  Km 
altitude.  Figure  2,  bottom  shows  ion  densities  at  finer  (1-s) 
temporal  resolution.  At  the  western  wall  of  the  depletion, 
ion  density  decreased  by  four  orders  of  magnitude  from 
5.104  to  ~5  cm-3  in  30  seconds  (~200  km  in  track).  The 
locations  of  the  dawn  terminator  at  the  altitudes  of  the 
satellite  (09:32  UT),  and  the  E-layer  (~09:35  UT)  are 
indicated.  C/NOFS  entered  the  depletion  shortly  before 
crossing  the  E-region  terminator.  Figure  2  also  shows  the 
presence  of  large-amplitude  irregularities  within  the  deple¬ 
tion.  While  the  west  wall  is  extremely  sharp,  the  east  side 
shows  a  gradual  upward  slope,  due  to  progressive  refilling 
of  the  ionosphere  from  below. 

[5]  The  vertical  component  of  the  V  =  (E  x  B)/B2 
plasma  drift  is  also  plotted  in  Figure  2.  The  drift  turned 
upwards  about  600  km  to  the  west  of  where  C/NOFS 
entered  the  depletion  and  remained  so  across  it.  Three  large 
(~200  m/s)  upward  spikes  appear.  The  first  occurred 
~1.3  minutes  before  the  ion  density  minimum.  The  third 
is  collocated  with  the  depletion’s  west  wall.  The  drifts 
presented  in  Figure  2  include  the  structured  drifts  associated 
with  the  dawn  plasma  cavity,  and  are  enhanced  compared  to 
the  background  drift.  During  this  period,  we  estimate  that 
the  peak  background  drift  was  ~80  m/s.  This  value  was 
successfully  used  to  model  the  background  density  deple¬ 
tion  by  Su  et  al.  [2009],  Although  we  do  not  understand 
well  the  cause  of  these  enhanced  drifts  and  their  relation  to 
the  density  cavity,  the  drifts  are  not  large  compared  to 
upward  drifts  of  classical  spread-F  depletions,  which  some¬ 
times  are  observed  with  speeds  of  500  m/s  or  even  1  km/s. 
Note  also  that  the  conductivities  during  this  period  are  so 
low  that  high  electric  fields  can  be  maintained.  There  is  no 
plasma  to  short  them  out.  In  addition,  the  reported  E-fields 
are  similar  to  the  enhanced  fields  observed  near  sunrise  and 
reported  by  Aggson  et  al.  [1995]  (see  event  C  of  Aggson  et 
al.  [1995,  Figure  1]). 

[6]  Ionosonde  measurements  on  this  day  from  Jicamarca 
corroborate  the  C/NOFS  VEFI  measurements.  The  low 
plasma  density  precluded  using  ISR  measurements,  but 
the  ionosonde  data  can  be  used  to  estimate  upward  drifts. 
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Figure  1.  Dawn  plasma  depletion,  orbit  915,  Jun  17,  2008. 
(top)  1-min  integrated  ionospheric  density  (black),  1-s 
integrated  density  (green),  satellite  altitude  (blue).  Middle 
panel:  FFT  from  the  512  Hz  PLP  data,  (bottom)  Map  with 
satellite  trajectory  (red)  and  magnetic  equator  (green). 
Thicker  dotted  lines  indicate  when  C/NOFS  was  in 
darkness. 


Between  09:45  and  10:00  UT  (~05  LT),  the  F-layer  peak 
altitude  increased.  The  velocity  inferred  from  this  rise  was 
~  1 5  m/s.  Then,  between  10:00  and  10:05,  the  upward 
velocity  abruptly  increased  to  125  m/s.  After  10:05,  it  was 
not  possible  to  accurately  measure  the  F  peak,  because  the 
density  was  too  low  to  allow  for  the  o-mode  return  to  be 
seen.  Using  he  x-mode  return,  we  infer  that  the  upward  drift 
may  have  remained  high,  possibly  reaching  50  m/s,  until 
10:35  UT.  The  ionosonde  return  is  an  average  over  the  large 
region  covered  by  the  antenna  beam,  necessarily  much 
lower  than  a  point  measurement  from  a  satellite.  Therefore 
the  ionosonde  data  confirm  the  presence  of  large  upward 
drifts  at  dawn  on  that  day  and  are  consistent  with  the  C/NOFS 
satellite  VEFI  in  situ  measurements. 

[7]  Figure  3  displays  ion  densities  sampled  during  a 
nearly  simultaneous  pass  of  the  Defense  Meteorological 
Satellite  Program  (DMSP)  satellite  F17  at  850  km  near  the 
dawn  meridian.  FI 7,  a  sun-synchronous  satellite,  crossed 
the  magnetic  equator  at  ~09:24  UT,  about  10  minutes 
before  C/NOFS’  first  encounter  with  the  dawn  depletion. 
The  plotted  measurements  show  that  FI 7  detected  a  latitu- 
dinally  broad  depletion  with  embedded  irregularities.  An 
examination  of  DMSP  passes  reveals  that  encounters  with 
similar  depletions  at  dawn  are  not  unusual  during  this 
period. 

[8]  C/NOFS  crossed  similar  depletions  at  all  altitudes 
between  perigee  and  apogee.  Figure  4  provides  a  second 
example  of  a  dawn  plasma  depletion.  It  occurred  at  23:09 


UT  while  the  satellite  was  at  650  km  in  the  Indian  longitude 
sector.  A  Fast  Fourier  Transform  (FFT)  of  ion  densities 
sampled  at  512  Hz,  represented  as  color  spectrograms  in  the 
middle  panel,  indicates  that  irregularities  appeared  within 
the  dawn  depletions  but  not  elsewhere  during  this  pass. 

[9]  This  orbit  was  selected  because  CHAMP,  a  polar- 
orbiting  satellite  at  325  km  altitude,  traversed  the  same 
dawn  depletion  ~20  min  earlier.  Figure  5a  displays  in-situ 
ion  densities  measured  by  CHAMP  during  ~  I  orbit  (red 
line).  International  Reference  Ionosphere  (IRI)  plasma  den¬ 
sities  along  the  CHAMP  trajectory  are  also  shown  (green 
line)  for  comparison.  CHAMP  crossed  the  magnetic  equator 
at  ~05  and  1 7  LT.  The  pronounced  density  minimum  at  the 
center  of  Figure  5a  coincides  with  the  satellite’s  crossing  of 
the  magnetic  equator  on  the  morning  side  at  22:47  UT.  The 
density  fell  from  105  to  2  x  103  cm-3.  Note  that  IRI  failed 
to  reproduce  this  depletion,  predicting  densities  well  above 
all  measured  values. 

[10]  Figure  5b  summarizes  characteristics  found  in  the 
present  example.  The  depletion  is  located  just  to  the  west  of 
the  E-region  terminator.  The  dimensions  of  this  dawn 
depletion  deduced  from  the  roughly  orthogonal  C/NOFS 
and  CHAMP  orbits  are  14°  in  longitude  and  50°  in  latitude. 
The  oval  sketched  in  this  map  approximates  the  size  of  the 
depleted  ionosphere,  and  illustrates  its  large  extent. 

[11]  Almost  identical  dawn  depletions  were  seen  by 
CHAMP  during  seven  consecutive  orbits  on  this  day.  Each 
time,  the  depletion  occurred  as  the  satellite  crossed  the 
magnetic  equator  and  very  close  to  the  E-region  terminator. 
Electron  temperatures  (not  shown)  measured  by  CHAMP 
within  these  depletions  were  1,000°,  as  predicted  by  IRI. 
Several  days  earlier,  when  the  satellite’s  orbit  was  near  the 
06  LT  meridian,  electron  temperatures  peaked  near  3100° 
during  the  equatorial  crossing. 

[12]  Ion  densities  measured  by  CHAMP  at  the  magnetic 
equator  were  averaged  in  one-hour  bins,  from  03  to  07  LT 
(L.  McNamara,  private  communication,  2009).  Morning 
side  density  depletions  are  deepest  in  the  05  LT  bin.  In 


Figure  2.  High  time-resolution  ion  density  and  plasma 
drift,  17  June  2008. 
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Figure  3.  DMSP  F17  ion  density,  17  June  2008. 


September  2008  they  were  deepest  in  two  longitude  sectors: 
30°  to  110°,  and  310°  to  350°. 

3.  Discussion 

[n]  To  recapitulate,  the  C/NOFS  satellite  repeatedly 
observed  deep  plasma  density  depletions  close  to  its  cross¬ 
ings  of  the  E-layer  terminator.  They  are  associated  with 
irregularities  and  were  observed  at  all  C/NOFS  altitudes 
(400  to  850  km).  Similar  depletions  were  crossed  by  the 
polar-orbiting  DMSP  and  CHAMP  spacecraft.  Closely 
coincident  passes  of  C/NOFS  and  DMSP/CHAMP  indicate 
that  the  depletions  cover  about  14°  in  longitude  and  50°  in 
latitude. 

[14]  Preliminary  results  show  that  occurrence  rates  of 
these  dawn  depletions  depend  on  season  and  longitude. 
For  example,  between  September  10  and  15,  2008,  dawn 
depletions  were  seen  in  77%  of  C/NOFS  passes.  However, 
significantly  fewer  events  were  observed  between  November  1 
and  6,  2008  when  dawn  depletions  were  recorded  in  only 
15%  of  the  orbits.  In  all  seasons,  depletions  appear  more 
frequently  and  are  deeper  in  the  America-Africa  and  India- 
Indonesia  longitude  sectors.  The  presented  data  were  ac¬ 
quired  during  magnetically  quiet  times.  The  Kp  index  was 
2+  and  1 — ,  for  the  June  and  September  examples,  respectively, 
and  remained  low  throughout  the  entire  two  days.  Conditions 
were  generally  quiet  from  10  to  15  Sept  with  Kp  <  3, 
except  for  the  last  9  hours  of  the  1 5th  when  Kp  =  4. 

[15]  Scherliess  and  Fejer  [1999]  indicate  that  upward 
plasma  drifts  appeared  in  their  recorded  satellite  and  ground 
data  and  were  subsequently  incorporated  into  their  clima¬ 
tological  model.  For  example,  scatter  plots  in  their  Figures  1 
and  2  show  many  instances  of  upward  drifts  in  the  dawn  LT 
sector.  Their  plots  show  peaks  of  20  and  40  m/s,  respec¬ 
tively.  These  represent  background  drifts  and  are  of  similar 
magnitude  as  the  C/NOFS  drifts  reported  here.  Interestingly, 
their  Figure  4  shows  large  discrepancies  between  the  model 
and  average  drift  values  at  night  after  ~20:00.  As  indicated 
above,  large  upward  drifts  at  dawn  were  also  deduced  from 
the  Jicamarca  ionosonde.  Thus,  the  C/NOFS  measurements 
are  not  unique  in  showing  significant  upward  drifts  around 
05  LT. 

[16]  Retterer  et  al.  [2005]  developed  the  Physics  Based 
Model  (PBMod)  to  assist  C/NOFS  analysis.  Su  et  al.  [2009] 
ran  PBMod  twice  for  the  entire  June  17day:  first,  using 
VEF1  measurements  and  then,  using  the  climatological 
electric  fields  of  Scherliess  and  Fejer  [1999]  as  drivers. 
As  mentioned  above,  the  PBMod  simulation  using  VEFI 


inputs  reproduced  the  observed  dawn  density  depletion.  It 
also  showed  the  height  of  maximum  density  HmF2  reaching 
~900  km  with  well  separated  Appleton  anomaly  peaks.  The 
climatological-driven  simulation  reproduced  neither  feature. 

[17]  This  paper  is  not  the  first  to  report  dawn  depletions. 
Burke  et  al.  [1979]  described  dawn-sector  depletions  ob¬ 
served  by  the  DMSP  F2  satellite  under  solar  minimum 
conditions.  However,  electric  fields  were  unavailable  at  that 
time.  Oya  et  al.  [1986]  inferred  plasma  density  depletions  at 
dawn  from  responses  of  an  impedance  probe  on  the  Hinotori 
satellite  during  the  1981-1982  solar  maximum.  ROCSAT-1 
also  crossed  density  reductions  near  sunrise  during  solar 
maximum  storms  [Su  et  al.,  2009],  Fejer  et  al.  [1999] 
mention  density  depletions  close  to  sunrise  and  attribute 
them  to  the  stormtime  disturbance-dynamo. 

[is]  We  are  now  left  with  the  challenge  of  explaining  the 
presence  of  an  eastward  E  field  at  dawn.  The  electric  field 
could  be  generated  by  a  dynamo  mechanism  driven  by 
winds  and/or  gravitational  drifts  [Eccles,  2004],  or  could  be 
due  to  over-shielding  by  region  2  Birkeland  currents.  The 
winds  that  drive  the  dynamo  could  be  due  to  geomagnetic- 
disturbance  dynamo  or  a  traveling  atmospheric  disturbance 
triggered  by  the  passage  of  the  terminator  [Fujiwara  and 
Miyoshi,  2006].  Non-migrating  tides  in  the  equatorial  ion¬ 
osphere  [Immel  et  al.,  2006]  and  the  strength  of  the 
geomagnetic  field  may  explain  part  of  the  longitudinal 
variation.  If  the  ultimate  cause  is  a  dynamo  produced  by  a 
zonal  flow,  then,  in  a  way  similar  to  the  pre-reversal 
enhancement  that  occurs  at  dusk,  we  can  invoke  the 
requirement  for  current  continuity  across  the  discontinuity 
of  conductivity  at  dawn  to  explain  the  generation  of  electric 
fields  in  the  orthogonal  direction,  leading  to  vertical  plasma 
drifts. 

[19]  Concerning  the  generation  of  these  fields  by  geo¬ 
magnetic  activity,  Burke  et  al.  [2009]  and  Su  et  al.  [2009] 
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Figure  4.  Same  as  Figure  1,  but  for  10  September  2008. 
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Figure  5.  (a)  CHAMP  data,  10  September  2008  with 
density  (red)  and  IRI  model  (green),  (b)  Sketch  of  the 
depleted  ion  density  region  in  the  dawn  depletion,  10 
September  2008.  Both  CHAMP  and  C/NOFS  see  dawn 
depletions  when  the  E  region  below  the  satellite  is  still  in 
darkness,  but  very  close  to  the  terminator.  Oval  illustrates 
the  depletion  size  ~50°  x  14°. 


argue  that  the  event  of  June  17  may  be  due  to  the  disturbance 
dynamo  or  over-shielding  by  Region  2  currents.  A  solar  wind 
disturbance  reached  the  Earth  on  June  14.  Although  it  did  not 
drive  a  magnetic  storm,  it  did  cause  irregular  increases  in  the 
AE  index  during  the  following  days.  AE  reached  ~700  nT  at 
05  UT  on  June  16.  The  authors  of  these  two  papers  regard  the 
Figure  1  depletion  as  a  remnant  of  a  wider  depletion  that  had 
formed  earlier.  Nevertheless,  this  remnant  cannot  be  consid¬ 
ered  a  dead  depletion.  It  was  kept  alive  and  nurtured  by 
continuing  large  upward  drifts.  It  is  clear,  however,  that  many 
dawn  depletions  observed  during  C/NOFS  passes  cannot  be 
regarded  as  remnants  of  prior  geomagnetic  activity.  They 
appear  spontaneously  as  isolated  events  that  start  at  about 
05  LT,  with  no  discernible  antecedent  signatures. 


4.  Conclusion 

[20]  We  have  shown  that  deep  pre-sunrise  depletions 
associated  with  upward  plasma  drifts  were  frequently  ob¬ 
served  by  C/NOFS.  Coincident  DMSP  and  CHAMP  data 


provide  an  estimate  of  the  extent  of  these  depletions.  Using 
the  electric  field  measured  by  VEFI,  PBMod  has  reproduced 
the  observed  features  of  the  depletions.  We  considered 
several  possible  reasons  for  the  generation  of  the  observed 
upward  drifts.  As  part  of  our  future  work  we  plan  to 
investigate  in  more  detail  possible  explanations  for  these 
depletions  at  05  LT.  In  addition,  as  C/NOFS  measurements 
accumulate,  we  will  derive  statistical  models  of  this  effect  in 
order  to  understand  how  and  why  the  dawn  depletions  and 
corresponding  E-fields  vary  with  longitude,  local  time, 
season,  magnetic  activity  and  solar  flux. 
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